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Abstract 



We study the phase transitions of a random copolymer chain with 

quenched disorder. We apply a replica variational approach based 

•^ ' on a Gaussian trial Hamiltonian in terms of the correlation functions 

of monomer Fourier coordinates. This allows us to study collapse, 
phase separation and freezing transitions within the same mean field 
^S| ' theory. The effective free energy of the system is derived analytically 

^^ I and analysed numerically. Such quantities as the radius of gyration or 

the average value of the overlap between different replicas are treated 

as observables and evaluated by introducing appropriate external fields 

'jrt ' to the Hamiltonian. We obtain the phase diagram and show that this 

Cj I system exhibits a scale dependent freezing transition. The correlations 

I . between replicas appear at different length scales as the temperature 

' O ' decreases. This indicates the existence of the topological frustration. 

a' 
o 

P. '. 1 Introduction 

>. 

^ , Various models of random copolymers have been extensively studied using 

H i theoretical [Q, §, ^, |, |5|, H and numerical |0, |, ^ methods. It has been 

found that even simple models exhibit a rich variety of conformational states 
leading to quite a complicated phase diagram. 

Most of the theoretical results have been obtained within the replica 



approach [^0[ for directed polymers in random media [12, |l^, 14 1 and for 



self-interacting copolymers under the assumption of the ground state domi- 
nance and constant density approximations for Gaussian disorder [§, ^, |^. 
Computer simulations [^, |^, ^] have qualitatively confirmed some of the the- 
oretical predictions. 

In this paper we study the complete phase diagram of a random copoly- 
mer chain with attractive and repulsive units. It is believed that such studies 
have application in the field of protein folding and related subjects |jll|. Thus 
we apply a Gaussian variational approach based on a generic quadratic trial 
Hamiltonian in terms of monomer Fourier coordinates of replicas, which is 
an alternative to the standard one using the density overlap function as 
the replica order parameter and constant density approximation. A similar 
approach was applied to directed polymers and manifolds in random media 
[|13| , P^ ]. The method has the advantage that it allows us to incorporate fluc- 
tuations of the density, determined self-consistently, and to study collapse, 
phase separation and glass transitions within the same mean field theory. 

Note that the self-consistent method has been applied to various other 
problems at equilibrium fig , 17 1 and in kinetics ||l^ of homo- and block 



copolymers and to directed polymers in random media [13, 14]. 

It is convenient to represent the Hamiltonian as the sum of two parts, 
H = Hhom + Hiiis- The first is purely homopolymeric and the second is 
related to the quenched disorder. Thus, 

PH,^^=u^dsm^y (1) 



ds 



+ E 7^^ / {d^o-dsm} n '^(K^o) - r(..)) 



1 /-W i-lV 

Hdis = ^ / dsds' u{\s, K') S{ris) - r{s')) , (2) 

^ JO Jo 

where the connectivity constant is k = 3/~^, with I being the bond length. 
Here the u„ denote the virial coefficients and the curly braces indicate that 
the terms with coinciding indices are excluded. The matrix of two-body 
interactions may be written |]l| via the quenched random variables {As}, 

'U2(As, Xs') =v + a{Xs + Xs') + xXsXs', (3) 

where the parameters v, a and x may be expressed in terms of, for example, 
UAA, ubb and uab in the case of only two monomer types. The average 



over the quenched disorder is performed by introducing replicas [10| 

PFguench = -logZ({AJ) = - hm n-1 • (^Z"({A,}) - l) , (4) 

Z^{{Xs})= J Drtexp{-(3Hn) (5) 

where Z, Hn are respectively the partition function of the system and the ef- 
fective Hamiltonian in the replica space. Here the bar stands for the average 
over disorder. 

For simplicity we assume the distribution of {As} to be Gaussian with 
mean value Aq and variance A^, 

mA,})=np^»p(-(^y£)!). (0) 

However other distributions, for example the binary distribution, can be 
considered. Averaging over the quenched disorder produces a compact ex- 
pression for the effective Hamiltonian in the replica space, 

r«(s,si)EEr'^(s)-r'^(,i), (8) 

where we assume x = 0, Aq = and that the Hamiltonian is regularized 
by introduction of a short distance cutoff. This choice of parameters in 
the matrix of interactions corresponds to the amphiphilic model, in which 
monomers of j4-type attract each other and monomers of -B-type repel each 
other. We note here that the case x 7^ requires some additional treatment, 
for example application of the perturbation expansion of the matrices M~^ 
and logM, where 

M{s, s') = 6{s - s') + PA\ Yl ^i^"(') - r"(^')) • (9) 



2 Variational approach 

It is useful to introduce the Fourier coordinates of the monomer positions. 
For ring polymers the appropriate Fourier transform is defined as 



1 [^ , f-i2Trsq\ , , , , ^ fi2Trsq\ 

— J dsexpl — — — lr(s), r(s) = tq + 2^ exp ( ^;^ Ir^ , (10) 



and the generalization to the case of open chain is straightforward. The 
standard Gibbs-Bogohubov variational principle gives the estimate for the 
free energy, 

F„ = Fo + ((F„-Fo))o. (11) 

We choose the trial Hamiltonian as a quadratic form, which is standard 
in the polymer literature [12, 13 1, but with g-dependence of all variational 
parameters 

PH, = \Y.ll^aM^%A^ (12) 

a,fe <? 

where the effective potential matrix is to be taken non-diagonal in the replica 
indices and diagonal in the chain index due to the assumption of transla- 
tional invariance along the chain after integrating out the disorder. We note 
here that formally one has to add a small mass term ^5ab to the effective 
potential V^^ in (02) for regularisation 113]. The limit /i ^ is assumed to 



be taken at the final stage of calculations. 

Averaging ^ over the statistical ensemble (12) we obtain the effective 
free energy of the following form 

q a q 



m- 

,2 



'^(m + l)!(27r)3W2^ 



2(2-) . „^, 

Here we introduced A = (3Aa and matrices of spatial D^^^^^ = i((r"(si) — 
r (s2))^)o and Fourier Tg = ^(r" rq)o monomer correlations 

The procedure of performing the limit n ^ is well known from studies 



of spin glass systems IT^, 12] . To do so we assume that !F^^ is a Paris! 



type hierarchical matrix parametrised by {fig,!Fq{u)}, u G [0,1], where Tq 
and J-q{u) are correlation functions for monomers from the same replica 
and from different replicas respectively. Finally, we obtain the effective free 



energy per replica in the limit n — > 0, which is convenient to represent as a 
sum of two terms, 



f5F,ff = lim(/?F„/n) = PN[ MJ',] + MT„T,{u)] ) . 



The first term, /q, coincides with the annealed effective free energy, 
PNfo 



-^ETrlog^. + ^E^^^. 



N 



+ t n:\~J'':'Z I {dso..ds,}x..r'^' 



fr{{l + l)!(27r)3'/2 
where Xi...^ are determinants of tti x ?7i-matrices 



(14) 



(15) 



The second term, /i, is determined purely by the quenched disorder, 



pNfi 



E 



Log (1 - (7^,)) + 



TZgiO) 



+ 



^Log 
u 

A2 



1 - (^5) 

l-(7^,)-[7^,](^) ^ 
1 - {n,) 



(16) 



2!(27r)6/2 7 ^^°^^^^'^V ^'"['^^o^i'^^i^^ -^sosisal"")] ^^^ 



where by definition for any arbitrary matrix a: 

"1 



a 







dua{u) , [a](tt) 



dv a{v) + ua{u) 



and 



0<7^g(^l) = ,Fg(n)/J^g<l, 



(17) 



(18) 



Here we introduce TZq{u), which gives the measure of correlations between 
different replicas relatively to correlations in the same replica and allows 
the division on "annealed" and "quenched" parts (|l4 )- The effective free 
energy has to be minimized with respect to Tq and maximized with respect 
to Tg{u) with the constraints < TZg{u) < 1. The latter physically means 
that monomer coordinates' correlations between different replicas are always 



weaker than those in the same rephca. Alternatively these constraints can be 
derived if we require the matrices Vab{q) and J^^^ to have positive eigenvalues 
for arbitrary integer n, n > 1, and all non-diagonal elements of the matrices 
Vab{q) to be negative to match the sign of the corresponding non-diagonal 
terms in the effective Hamiltonian. 

Usually one is interested in such quantities as radius of gyration Rg, end- 
to-end distance, phase separation order parameters or the average value Q 
of an overlap between different replicas Qab^ for example, 



R^ = l/2iV-2 f dsds'{{r{s) - r(s')) 



2^0 



Qab = l^N-^ J ds{6{r-{s) - r^{s)))o , (19) 

6Rl = l/2N-^ J dsds'{Xs{r{s)-r{s'))^)o . 

All these can be treated as observables and evaluated by introducing ap- 
propriate external fields to the partition function. The model we consider 
exhibits a strong tendency for phase-separation. Therefore, the parameter 
^ = 6R^/R^ can be considered as a measure of phase separation, where 6R^ 
is the difference of the square radii of gyration in the case of two letter copoly- 
mers with equal concentrations of both monomer types, 6Rg = R^j^ — R^b- 
Introducing the term i^iSR^ into the Hamiltonian one obtains 



li,n- 



■*o d^ dn \ 



8Rl = - hm^— — (Z"(/.)> (20) 



^ ^ / ^„_^„,^„_ ■50S1S2 X-^SOSI 



J„ J„ J„ ^gOSlS2 \^S0S\S2l 

4 (2vr)3/2Ar2 7 ^^^od^id.s ^^^y 

To determine the freezing transition one may study the change of the 
average overlap among replicas Q as well as the eigenvalues of the matrix 
of the second derivatives of the free energy (|l^ with respect to all Tq and 
Tq{u), Hessian T. The solution we are looking for is a saddle point with the 
positive eigenvalues of T corresponding to J^q and negative ones correspond- 
ing to J-q{u). As one can see from the numerical analysis, the solution with 
all J~q{u) = 0, no correlations among replicas, is not a saddle point of the 
type mentioned above for the variance of the disorder A larger than some 
particular critical value Aj . The latter can vary with a change in the values 
of the virial coefficients. 



3 Numerical Analysis 

To analyze system numerically one has to substitute integration by summa- 
tion in dl^), see Appendix A for details. In principle we can evaluate the 
effective free energy for any finite number of steps in the Parisi scheme. In 
practice this means the introduction of a large number of variational param- 
eters. If one has a polymer of length A'^ monomers and applies the fc-step 
Parisi scheme, then one has to analyze a set of approximately N{k + 1) non- 
linear algebraic equations. Therefore, we apply a minimal one-step Parisi 
scheme, which still incorporates all the main features of the exact solution 
and allows one to locate the onset of the transition to phases with broken 
replica symmetry if such exists. Thus we approximate J-q{u) by the step 
function 



Tq{u) 






(0) 



q u < X, correlations between replicas from different groups, , , 

(1) \^^) 

J-q u > X, correlations between replicas from the same group. 



in one-step approximation: 



/3^/i = -|Ef-i°g(i-(i-^)^^ 



(1) 



E 



x-l 



+ 



q L 

A2 



^o,{i-ni^^)^ 



x7^w 



n. 



(0) 



2!(27r)6/2 



E 

{tjk} 



(1 - x)7^, 



(1) 



xTi, 



(0) 



(22) 



q!;1 = aA. - ^f , ^gi = E 4^^^^ ' ^ = 0' 1 



(23) 



In the following we restrict ourselves to the choice of parameters: fc^T = 
1, / = 1, which provide the scaling units of energy and length; and we set 
Um = ioT m > 5. Note that the zero-modes ^o and T^q ' do not contribute 
to the mean energy of the system describing diffusion of the center of mass 
of the polymer chain. 

The system exhibits two different types of behavior above and below the 
freezing transition. Above the freezing temperature, the monomer coordi- 
nates correlation function for monomers from the same replica, jFg, closely 
coincides with that for annealed disorder, whilst the monomer coordinates 



correlation functions for monomer from different replicas, TZq and TZq , 
possess only trivial solutions, namely: TZq = TZq = 0. The coil and globu- 
lar states are similar to the homopolymer coil and liquid-like globule states 
respectively. The trial Hamiltonian is diagonal in the replica space and the 
free energy differs from the annealed one by the term 

{momim2} 

which remains regular in the whole range of the parameter space and be- 
comes negligibly small in the region of small A. 

The situation changes around some critical value of the dispersion of 
disorder Aj, at which point the solution yields nontrivial values for the 

correlation functions TZq and TZq and the parameter x. Thus, the onset of 
the freezing transition can be seen clearly from the behavior of both x and 
the set of TZq ' , where they start to change rapidly. The typical behavior 
of the correlation functions for monomers from different replicas can be seen 
in Figs 1 and 2. 

Fig. 1 presents the correlation function for monomers from the same 
group of replicas, 7^i=i, and from different groups of replicas, 7^1=1, and 
the parameter x versus the dispersion of disorder A at fixed values of the 
other parameters in the free energy (|l^). Strictly speaking the solution 
below freezing is valid only in the vicinity of the transition line, though we 
draw it for a wide range of the dispersion of disorder, keeping in mind that 
more steps in the Parisi scheme are required to describe the properties of the 
frozen phase. After the freezing transition the coordinates of monomers from 
different replicas become correlated which results in a decrease of the entropy 
of the system. The symmetry among replicas is spontaneously broken, the 
correlations between the positions of the monomers from the same group of 
replicas are higher than those between the positions of the monomers from 
different groups. 

Fig. 2 presents the profile of the replica's correlation functions TZq and 
TZq in the chain index q. There are no correlations for large q, q < N/2 
for ring polymers, between different replicas in the vicinity of the freezing 
transition although these become non-zero as A increases. To locate the 
onset of the freezing transition only correlations on the first mode, TZ\ 
and 7^]^ , and the parameter x are needed. This clearly indicates the scale 
dependence of the freezing due to the topological frustration in the system 



and is in agreement with the main result of Jl8| . 

The resulting phase diagram of the random copolymer is represented in 
Fig. 3. It shows that the collapse can occur in two different ways. For small 
values of the variance of the two body interaction, A, the coil collapses to 
the liquid-like globule as the effective two body interaction changes from 
repulsion to attraction and U2 becomes negative. For large A the effective 
three body interaction, {u^)(,ff = u^ — 3A^, changes sign from positive to 
negative and collapse occurs due to the three body effects [^ , even when -112 
remains positive, but small. There is also a transition from the liquid-like 
globule to the frozen globule (line (F) in Fig. 2). For very large A there 
is an indication that the copolymer chain collapses directly to the frozen 
phase and that the phase diagram may exhibit a multi-critical point where 
all three phases, coil, liquid-like and frozen globule coexist. 

The freezing transition line on the Fig 3 scales approximately as Aj ~ 
1^2! ''^ for large negative tt2, where 7 ~ 0.66 it 0.03. The freezing transition 
line is nearly independent of the third virial coefficient in the region where 
^3 <C n4 and increases linearly in ^3 when ti3 > ^4. In a more dense globule, 
i.e. given a bigger negative U2 or smaller u^, chains undergo a transition from 
a liquid-like to a frozen globule at smaller values of A. Both the collapse 
transition and the freezing transition lines are essentially independent of the 
size of the chain for sufficiently large degree of polymerization A^. 

In Fig. 4 we present the phase separation parameter 6R^ versus the 
dispersion A, 

. _ 3 A-A ^ ^,-xVf^-{l-x)V^J^ 



9 A (0^\i/2 AT2 Z^ 



4 (27r)3/2iV2 ^ ^5/2 

^ ' ^{3k} ^jk 

This parameter increases linearly in A in the liquid-like globule phase, but 
remains constant in the frozen phase, whilst ^ reaches a maximum at the 
freezing transition and starts to decrease slowly in the frozen phase due to 
the increase in the radius of gyration Rg . This means that freezing prevents 
the system from further phase separation and therefore, for example, the 
observable 6R'l and its first derivative with respect to A can also be used to 
locate the onset of the freezing transition. 



4 Discussion and Conclusions 

In this paper we have apphed a Gaussian variational approach to study 
the phase transitions of random copolymers. The essential feature of this 
approach is that it is based on a trial Hamiltonian in terms of monomer 
coordinates rather than density variables. For this reason it may be used 
to describe a polymer with arbitrary Gaussian disorder across the whole 
range of the collapse transition. Our equations reduce to the self-consistency 
equations for a homopolymer at equilibrium as the variance of the second 
virial coefficient vanishes. The most important result is that the freezing 



transition in the system is scale dependent |18|. 

We have performed a numerical analysis of the free energy in the one-step 
Parisi scheme. We note here that, firstly, more steps in the Parisi scheme 
may be required to describe the property of the frozen phase and, secondly, 
the replica symmetrical solution brings qualitatively the same results. 

It is believed that the model we have considered can provide a basis 
for understanding conformational states of biomolecules such as proteins 



[11|. The results for the freezing transition and the phase separation are in 



agreement with those obtained for this model by other methods [23]. We 
note that the results we have obtained are valid also for binary distribution 
of disorder, provided that higher virial coefficients remain positive. However, 
the whole diagram in this case will incorporate new features, e.g. the explicit 
dependence of the freezing transition on the chemical composition. 

We also note that one of the important drawbacks of the theory is that 
the Gaussian method yields an incorrect Flory exponent for a chain in a 
good solvent. However, this problem can be resolved in different ways and 



we refer the reader to the discussion of this point in |22]. 

Finally we may comment that despite the limitations, the current ap- 
proach offers one of the few consistent approaches to the global "phase" 
diagrams of random amphiphilic copolymer chain |^, p^ ]. 

The authors thank A. Gorelov, E. Timoshenko, I. Moskalenko and D. 
Hegarty for comments. One of the authors (A.M.) thanks the organizers of 
PeH 96, where part of this work was presented. 
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5 Appendix A 



To perform a numerical analysis we have to substitute integration by sum- 
mation, 



pN ^-1 a 2 N-1 

J d{siS2...Sm} "-^ Yl '— ^H ^^-^e ^^'^^ H •^^^^'^^T^^^'^) 



where {iii2---im} indecates that indeces 11,12, ■■■,im are not allowed to coin- 
cide. Thus we can rewrite ( [T5| ) and (p^) 

PNfo = -^ E log -^9 + ^'^N Y, T, sin2 ^ 
g q 

+ ^ (n + l)!(27r)3W2 Z- /i-« ' ^^'^^ 

n=l V ' / V / {mo...mn} 



2 g L 1 - {Ti-q) Jo w' 1 - (vc, 

A^ 

2!(27r)6/2 



{ jj/c} 



where A = /3aA and Xi..,m are determinants of m x m-matrices 
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Figure Captions 

Fig. 1 

Plot of the correlation functions TZ\ ' and the replica symmetry 
breaking parameter x vs the dispersion of disorder A for polymer 
with the degree of polymerization N = 50. Here U2 l~^ = —50, 
^3 l~^ = 60, n4 l~^ = 48. The dispersion Aj denotes the point 
of the freezing transition. 



Plot of the correlation functions 7^^ (left-hand side of the pic- 



Fig. 2 

lot of the 

ture) and TZq (right-hand side) vs the chain index q for polymer 
with the degree of polymerization N = 50. Here U2 /"^ = —50, 
^3 l~^ = 60, n4 1~'^ = 48. Lines (a) - (c) correspond respectively 
to the values of the dispersion of disorder A = 2 (near the freez- 
ing transition line), 3 and 4. For convenience we have drawn 
the plots only on one of the sides, and the correlation functions 

may be extended to the other one by the symmetry property: 
^(0,1) _^{o,i) 
/■Cg — /-CTv-g- 

Fig. 3 

The phase diagram of a random copolymer in variables of the 
dispersion. A, and mean value, U2 /~^, of the second virial co- 
eficient. Lines (C) and (F) correspond to collapse and freez- 
ing transitions respectively. This diagram has been obtained 
from the data for a polymer with the degree of polymerization 
A^ = 20, the third and the fourth virial coefficients are equal to 
M3 l~^ = 60, n4 /~^ = 48. No substantial qualitative changes are 
observed for larger A'^. Due to finite width of the collapse transi- 
tion for finite-size systems the line (C) is drawn by the maximal 
descent of the squared radius of gyration. 
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Fig. 4 



Plots of the phase separation parameter, 6Rg l^"^, vs the disper- 
sion of the second virial coeficient, A, for polymer with A^ = 20, 
U3 l~^ = 60, U4 l~^ = 48, Ao = 0, A = 5 and the values of second 
virial coeficient U2l~^ = —10 (diamonds), —40 (crosses), —70 
(quadrangles) and —100 (triangles). 
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